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Abstract 

We consider the semiclassical (zero-dispersion) limit of the one-dimensional focusing Nonlinear Schrodinger 
equation (NLS) with decaying potentials. If a potential is a simple rapidly oscillating wave (the period 
has the order of the semiclassical parameter e) with modulated amplitude and phase, the space-time plane 
subdivides into regions of qualitatively different behavior, with the boundary between them consisting 
typically of collection of piecewise smooth arcs (breaking curve (s)). In the first region the evolution of 
the potential is ruled by modulation equations (Whitham equations), but for every value of the space 
variable x there is a moment of transition (breaking), where the solution develops fast, quasi-periodic 
behavior, i.e., the amplitude becomes also fastly oscillating at scales of order e. The very first point of 
such transition is called the point of gradient catastrophe. We study the detailed asymptotic behavior 
of the left and right edges of the interface between these two regions at any time after the gradient 
catastrophe. The main finding is that the first oscillations in the amplitude are of nonzero asymptotic 
size even as e tends to zero, and they display two separate natural scales; of order 0{e) in the parallel 
direction to the breaking curve in the (x, t)-plane, and of order Obelus) in a transversal direction. The 
study is based upon the inverse-scattering method and the nonlinear steepest descent method. 
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1 Introduction and main results 



The focusing Nonlinear Schrodinger (NLS) equation, 

iedtq+]^e'dlq+\qW^Q, (1-1) 

where a; G M and < > are space-time variables, is a basic model for self-focusing and self-modulation e.g. 
it governs nonlinear transmission in optical fibers; it can also be derived as a modulation equation for 
general nonlinear systems. It was first integrated by Zakharov and Shabat [I] who produced a Lax pair for 
it and used the inverse scattering procedure to describe general decaying solutions (lim|2,|^o q(x^ 0) — 0) 
in terms of radiation and solitons. Throughout this work, we will use the abbreviation NLS to mean 
"focusing Nonlinear Schrodinger equation" . 

Our interest in the semiclassical limit (e — > 0) of NLS stems largely from its modulationally unstable 
behavior, as shown by Forest and Lee [2], the modulation system for NLS can be expressed as a set of 
nonlinear PDE with complex characteristics; thus, the system is ill posed as an initial value problem with 
Cauchy data. As a result, initial data 

g(x,0,e) = A(a;)e'*(")/^ (1-2) 

i.e. a plane wave with amplitude modulated by A(x) and phase modulated by ^(x) are expected to break 
immediately into some other, presumably more disordered, wave form when the functions A(x) and <i>(x) 
possess no special properties. 

However, in the case of an analytic initial data, the NLS evolution displays some ordered structure 
instead of the disorder suggested by the modulational instability, as shown by Fig. [l] see [3], [1] and [5]- 

Fig. [l] from [5j clearly identifies regions where different types of behavior of the solution q{x,t,s) 
appear. These regions are separated by some independent of e curves in the x,t plane that are called 
breaking curves or nonlinear caustics. Within each region, the strong asymptotics of q{x, t, e) can be 
expressed in terms of Riemann Thet a- functions to within an error term of order 0(e) that is uniform on 
compact subsets of the region ([B] for the pure soliton case and [7] for the general case). In this context, 
regions of different asymptotic behavior of q{x, t, s) corresponds to the different genera of the hyperelliptic 
Riemann surface whose Theta-functions enter in the asymptotic description. No error estimates for the 
asymptotics on the breaking curve were studied so far, although it was expected that the accuracy should 
be of the order 0{e^) at any regular point on the breaking curve. A visual inspection of Fig. [l] suggests 
that the behavior of q at the tip of the breaking curve is different from the behavior elsewhere on the 
breaking curve. 

The tip-point is called a point of gradient catastrophe, or elliptic umbilical singularity ([SJ).) 

The main goal of this paper is to analyze the leading order asymptotic behavior of the solution q{x, t, e) 
on and around the breaking curve except a vicinity of the gradient catastrophe point, and to obtain the 
corresponding error estimate. More precisely, we will examine a strip around the breaking curve between 
the genus zero and genus two regions, which has the width of order ©(elne). 
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1.1 Description of results 

The main goal of the paper is to provide a detailed asymptotic analysis in the zero-dispersion (or semi- 
classical) limit (e 0) in the neighborhood of the breaking curve, that is, the frontier in {x, i)-space that 
separates the region of smooth behavior {genus zero region) from the region of fast oscillations {genus 
two region), which is clearly seen on Fig. [l] 

The main highlight is an "universal" expression for the behavior of the first oscillations as we egress 
from the genus zero region into the genus two one. Here "universality" means that the expression does 
not depend upon the details of the initial data, or rather, it depends on it only through a few parameters 
that are explicitly computable. In order to give a taste of the expression in consideration we need to 
present a minimum of background notation, with more details contained in the body of the paper. 



The zero-dispersion limit of (1-1 1 can be addressed in terms of modulation equations involving one 
pair of complex conjugate Riemann invariants, say a, a. The invariant a — a{x,t) (which we choose in 
the upper half plane) depends on time/space according to the modulation equations (also called Whitham 
equations). The asymptotic behavior for the solution of the NLS equation in the genus zero region is 
then provided by |7| 

g(a;,i) - -6(a;,i)ei*(^'*)(l + o(l)), where (1-3) 

a{x, t) = a{x, t) + ib{x, t) 
d^<^{x,t) = -2a{x,t) 

and o(l) denotes some infinitesimal quantity in e. However, as we leave the genus-zero region, such 
approximation breaks down: the term that was infinitesimal o(l) becomes suddenly of finite order and 
a different approximation scheme must be employed. The modulation equations alone cannot detect 
such sudden change of behavior, since they are equations that are obtained from a formal manipulation. 
The modern approach to the problem requires the introduction of a special function, the (7-function; 
this is a locally analytic function in the complex spectral plane of the associated linear problem. It is 
defined in terms of the scattering transform of the initial data: without entering the details we may say 
that in our problem the initial data is encoded in a Schwarz-symmetrical analytic function fQ{z) (3/0(2;) 
generically has a jump discontinuity along the real axis). Taking into account the time (t) evolution and 
normalization with respect to x, one obtains 

f{z) := f{z- x, t) = fo{z) -xz- 2tz^ . (1-4) 

The g-function then is obtained as solution of a scalar Riemann-Hilbert problem (RHP) on a free 
boundary with the jump f{z) 

■^This means that the contour where the jump is has to be determined impUcitly along the way of the construction of 
the solution itself. 
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For the reader unfamiliar with the method we may compare the determination of the free-boundary 
to the determination of the steepest descent contour for Laplace-type integral. The main requirements 
(but this is a severely shortened list, see Section [2] for the details) is that 

9+{z) + 9-{z) ^ f{z;x,t) , ze-f^ (1-5) 

where 7™ is an arc of the said free boundary and the sign in the subscript denotes the left/right boundary 
value; moreover g is bounded at infinity and with some growth requirements at the endpoints (if any) of 
7m. The function g(z) = g(z; x, t), defined through the above mentioned RHP, must be supplemented by 
several inequalities that involve 

h{z; X, t) 2g{z; x, t) ~ f{z; x, t) , (1-6) 

and the sign distribution of its imaginary part 3(/i). All the functions /, g, h are Schwartz-symmetrical, 
in particular h(z;x,t) — h{z;x,t); this means that the sign-distribution requirements in the lower half- 
plane are the opposite of those in the upper half, C+. 

One of the main requirements is that S(/i) in C+ must be positive along the complement 7c of 
7m within the above-mentioned free-boundary, and must be negative on both sides of 7^, the latter 
condition implying that 3/i = on 7^ . What happens is that one of these particular requirements fails 
in certain regions of the (x, t) plane; this is where the zero-dispersion asymptotic behavior of the solution 
changes its behavior. The typical phase portrait is depicted in Fig. |4] The boundary between the two 
regions consists of points (x, t), for which 3(/i(z; x,t) = and simultaneously h'(z; x,t) = at some point 
z of the spectral plane; we denote by ri(x, t) the value of 2: € C-|- satisfying h' {z; x,t) — Q for all x, t in a 
vicinity of the breaking curve. Because of the analyticity of h this point 77 must be a saddle-point for 3/i 
and h can be written as 

h{z- X, t) = 5(x, t) - ^%^(z - r^Y + 0((z - 77)') ■ (1-7) 

(the factors are just for convenience). Note that the breaking curve is now given by 35(x,t) = 0. 
Generically, the second derivative coefficient C(x, t) does not vanish: this corresponds to the smooth part 
of the boundary between the two regions in the sketch of Fig. [4] The angular point of the boundary 
(where C(xo,io) = 0) is the so-called point of umbilical gradient catastrophe and will be dealt with 
in a forthcoming publication. 

The critical value S{x,t) G C plays the first major role in our story: indeed we have 

• The function S{x,t) is a local diffeomorphism for x G M, < e M+ with values in C ^ except at 



the point of gradient catastrophe (Lemma 3.2 1. 



The function 5(x, t) maps a tubular neighborhood of the breaking curve not containing the point 



of gradient catastrophe diffeomorphically on a tubular neighborhood of the real iS-axis (Thm. 3.1 ) 
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This means that we can use 55* as a parametrization in a transversal direction to the breaking curve, 
while using as a parametrization in the parallel direction. The paper focuses on the analysis of such 
a tubular neighborhood with a width of order ©(elne); more precisely we introduce the exploration 



parameter >c{x,t) (Def. 3.2 1 by 



S{x,t) 



(1-8) 



2 2 

The parameter >f is positive on the left breaking curve (negative on the right one, respectively) as we 
enter the oscillatory zone (in wavy lines in Fig. |4]); a finite value of x corresponds to an infinitesimal 
exploration. 




Figure 2: The relative amplitude of the corrected solution to the genus-0 solution. The plot is in the 
variables 9, k and the scale is not uniform. The parameters C, 77, etc. were obtained from an explicit genus- 
zero solution and are: a = 0.8120825800 + 1.251939783 i , ^ = -.5422408496 + 1.731859687 i , C = 
.2387253889 - .4667042683 i , Aq = -2.358599870 -I- 2.555370658 i , e = 10"!° , = 0.1 , to = 
0.1570482549. It is actually necessary to set such a small epsilon because for x e N we have a discontinuity 
of order -y^, to be expected since this is the order of our approximation and at the integers the correction 
we have is of the exact same order of the error. Note that since k measures distances in the scale 
ej lne| while in the scale e, the distance between consecutive "ranges", should be much longer than the 
transversal separation of the peaks. Also the successive ranges are relatively very thin (of order 1/| lne| 
relative to the longitudinal separation) . Therefore the typical size of the hills is e in all directions while 
the separation in the longitudinal direction is e\ lne|. This picture is consistent with a full-blown genus- 
two regime, where the solution is quasi-periodic at the scale e: as we progress into the genus-2 region the 
separation reduces to the natural scale e. If we were to plot this in the {x,t) plane, the only difference 
would be a linear change of coordinates and the "ranges" would be parallel to the breaking curve. Note 
that the picture is strictly periodic in the direction but only quasi-periodic in the x direction although 
it is almost not discernible in the picture: in fact the subsequent ranges are actually slightly shifted. 



The main finding of the paper is the fine structure of the "shock waves" along the breaking curve in 
the a;,t-plane, clearly visible on Fig. [T] They turn out to be not "uniform", as one might erroneously 
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conclude looking at Fig. [T] As shown on Fig. |2] each wave consists of the sequence of "peaks" or 
"teeth" interlaced by "holes" or "depressions" . It is remarkable that the mathematical description of 
this phenomenon is quite similar to the description of asymptotic of "colonization" of a new interval by 
zeroes of certain orthogonal polynomials in [SJE]. More precisely: 



• the genus zero approximate solution (1-3 1, continued into the genus-two region, gets a finite order 
correction (w.r.t. e) only for (positive) half— integer values of >r, which correspond to the rows 
of "teeth"; 



in between half-integers, the correction is of order O ^£2('i'>5*(^'^) i)^ ; 



• the correction term exhibits fast oscillations only in the direction parallel to the breaking curve, 
i.e., it depends on e'^^ and it is periodic in 

• all the remaining dependence of the correction is slow, namely, it depends on {x,t) at a 0(1) scale; 

• the formula for the correction term is valid with an error term of order 0{^); 

• the only dependence on e in the correction terms is through the expressions e^' ' and £^^^~^\ 
where K = max(0, \_>c\) with [>fj denoting the integer part of >f. 

The picture getting painted is that the first oscillation is of finite amplitude and the oscillations form 
"mountainous ranges" with fiatlands in between of order elne and peaks separated longitudinally by 
distances of order e. Possibly the most descriptive explanation is contained in Fig. [2] plotted over the 
(>r, t?) plane which, we remind, is diffeomorphic to a piece of the {x,t) plane around the breaking curve, 
with >c growing in the transversal direction (>f and are not necessarily orthogonal; in fact one could 



compute their relative angle easily from (3-3l) 

While the correction is strictly periodic in it is not periodic in xr, despite a superficial inspection 
of the picture: neither the formula for the correction nor the picture (at a careful inspection) are periodic 
in K. 

The actual formula for the correction term is significantly complicated and it is inconvenient to write 
it here: it is contained in Thm. |7.1| and -in quite explicit terms- in Thm. |7.2[ The main point that we 
observe at this stage is its universality, namely, the fact that it does not depend on the details of the 
initial data. To be more specific, it is an expression that contains only C, i?, >c, a, rj. 

Organization of the paper 

We start with Section [2] where the essentials about the inverse scattering transform are recalled following 
[7]. In this section we summarily go through the transformation of the associated Riemann-Hilbert 
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problem along the lines of the Deift-Zhou steepest descent method [TT] by introducing the g-function 
(Sect. \2.1\ and the transformation of the RHP (Sect. |2.2[ ). 

Section |3] we analyze the differential geometry of the breaking curve and establish the result that 



S{x,t) is a local diffeomorphism onto the spectral plane (Lemma 3.2 1. We also analyze the large-time 



behavior of S{x,t) along such breaking curve (Sect. 3.1) 



In Section |4] we start the detailed construction of the approximation to the RHP: the description 
is made in two steps, mainly for pedagogical reasons. In the first step we provide a leading order 
approximation to the problem and show that the error term fails necessarily to tend to zero as e — > for 



positive half-integer values of >€ (eq. 4-39 1 . This prompts the necessity of a higher order approximation 
scheme, contained in Sect. [5] 

The analysis is done firs on the right branch of the breaking curve and hence should be repeated for 
the left branch; this is done in Sect. [6j The finding is that the situation is different only in some setup, 
but the essence of the formulas and of the phenomenon is identical. 

The paper culminates in Sect. [7] where we provide the final formulae for the leading and subleading 
terms asymptotic behavior of the NLS solution (Thm. |7.1[ Thm. 7.2 1. We also provide a short discussion 
about the large-time behavior (along the breaking curve) of the solution. 



2 A short review of the zero dispersion Umit of the inverse scat- 
tering transform 



Given an initial data (potential) q{x,0) for the (1-1) that is decaying as x ^ ±oo, the direct scattering 



transform for the NLS ([1 ) produces the scattering data, namely: the reflection coefficient ro{z,e) and 
the points of discrete spectrum together with their norming constants (solitons). The time evolution of 
the scattering data is simple and well-known. Thus, to find the evolution of a given potential at some 
time t, one needs to solve the inverse scattering problem at the time t. Equivalently, we can stipulate 
that the initial data is assigned directly through the scattering data; thus, one can produce a solution 
to the NLS (1-1 1 by choosing some scattering data and solving the inverse scattering problem for i = 



(initial data) and for t > (evolution of the initial data) . The latter approach allows one to avoid solving 
the direct scattering problem and addressing many delicate issues associated with it (see [H] for more 
details) . Since we are interested in studying the generic structure of solutions of the NLS in the transition 



regions, it will be convenient for us to define a solution to the NLS (1-1 1 by its scattering, not by its 
initial, data. 



In considering the semiclassical limit of ( 1-1 1, one has to consider the semiclassical limit of the corre- 



sponding scattering transform. For the case of decaying potentials of type ( 1-2 ) , this hmit was discussed 
in |13j (inverse scattering) and |12j (direct scattering) , where it was shown to be a correspondence between 

1 



a{x, t) 



-<^^{x) + iA{x,t) 



(2-1) 
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on the potential side and 

/(z) = /o(z) ^xz^ 2tz^ (2-2) 

on the scattering side. Here f > is fixed and fo{z) has the meaning of the "scaled" logarithm of the 
reflection coefficient rQ{z,e) that corresponds to the initial data (1-2), that is, 



- liniero(2;,£) . 



(2-3) 



For example, the limiting fo{z) that corresponds to the potential cosh(a;)e^e' tanh^+iTr^ given by (see [7j) 



/o(z) = (l-z) [z-+ln(l-z) 



z In ; 



TT 

ln2+-e, 



when 3z > 0. 



(2-4) 



Moreover, if x{a) denotes the inverse function to a{x,t), at t fixed, then x(a) and /o(-z) are connected 
through some Abel type linear integral transform, see details in |12j . 

Since we have taken the perspective that we are starting with the scattering data of the form ro{z, e) := 
e-^fo{z) g]^all assume that ([T3]): 

• fo{z) is analytic in the upper half-plane and has a continuous limit on M; 

• fo{z) = 0{z) as z ^ oo; 

• there exists an interval (^_,^+) C M such that '^fo{z) > for z e (^_,//+) and 3/o(z) < for 
z e (-cx),/i_) U (m+,oo); 

• Qfo{z), z G M, has simple zeroes at fi± and its values are separated from zero outside neighborhoods 
of /i±. 

Due to the Schwarz-symmetry of the scattering data for the NLS, /o(z) is also a Schwarz-symmetrical 
function, i.e., /o(z) = fo{z); generically, Q/q^z) has a jump across the real axis. The above assumptions 
are not overly restrictive: it was shown that in the case of a solitonless potential the formula (2-3 1 yields 
such an /o(z) after some process of analytic extension, or "folding" of certain logarithmic branch-cuts 
onto the real axis. In fact, for our goals it is sufficient to replace the analyticity requirement by a weaker 
condition: the contours of the RHP for the ^-function (which will be introduced in the next section) 
always lie within the domain of analyticity of fo{z) in C-)-. Generically, ro(z) has a jump discontinuity 
along the real axis due to the discontinuity in Q/q. 



At any time t, the inverse scattering problem for (1-1) (in the solitonless case) with a fixed (not 
infinitesimal) e is reducible to the following matrix RHP. 

Problem 2.1 Find a matrix T{z) analytic m C \ M such that 



r+(z) = r_(z) 



roe 



rol^-Hl roe 

T{z) = l+-Ti + 0{z-^) , 



z e 



(2-5) 
(2-6) 
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(In the case with solitons, there are additional jumps across small circles surrounding the points of discrete 
spectrum.) Then 

q{x,t,s):^-2(T,)^^ (2-7) 



is the solution of the initial value problem (1-1) for the NLS equation 



The jump matrix for the RHP admits the factorization 



+ 1 r ' 




■ 1 r ' 




■ 1 " 


r 1 




1 




r 1 



(2-8) 



where r = r{z;x,t) = ro(2:)e (^*^^+^^) . 

Inspection of the RHP shows that the matrix T{z)* (where * stands for the complex-conjugated, 
transposed matrix) solves the same RHP with the jump matrix M(z) replaced by M~^{z) and hence 



Proposition 2.1 The solution T{z) of the RHP for NLS has the symmetry 

T{z){mr ^ 1 . 



(2-9) 



In order to study the dispersionless limit £ 0, the RHP ( 2-5 1-( 2-6 ) undergoes a sequence of transfor 



mations (that are briefly recalled in Section 2.2 1 along the lines of the nonlinear steepest descent method 
[HI [7] , which reduce it to an RHP that allows for an approximation by the so-called model RHP. The 
latter RHP has piece-wise constant jump matrices (parametrically dependent on x, t) and, in general, can 
be solved explicitly in terms of the Riemann Thcta functions, or, in simple cases, in terms of algebraic 
functions. The g-function, defined below, is the key element of such a reduction. 

2.1 The (yf-function 

Given fo{z), we introduce the g-function g(z) = g(z;x,t) as the solution to the following scalar RHP: 

1. g{z) is analytic (in z) in C \ 7m (including analyticity at oo); 

2. g{z) satisfies the jump condition 

g+ + g- = fa - xz - 2tz^ on 7„, (2-10) 

for a; e M and t > 0, and; 

3. g{z) has the endpoint behavior 

g{z) — 0{z — a) 2 -I- an analytic function in a vicinity of a. (2-11) 

Here: 



10 



• 7r?i is a bounded Schwarz-symmetrical contour (called the main arc) with the endpoints a, a, ori- 
ented from a to a and intersecting M only at 

• g± denote the values of g on the positive (left) and negative (right) sides of 7^; 

• the function /g = fo{z), representing the initial scattering data, is Schwarz-symmetrical and Holder- 
continuous on 7m. 

Taking into the account Schwarz symmetry, it is clear that behavior of g(z) at both endpoints a and a 
should be the same. 



Assuming /o and 7™ are known, the solution g to the scalar RHP (2-10 1 without the endpoint condition 

/(C) 



(2-11) can be obtained by the Plemelj formula 

R{z) 



2-111 



-dC , 



(2-12) 



where R{z) — ^y{z — a){z — a). We fix the branch of R by requiring that lim^^oo = 1. If foiz) is 
analytic in some region S that contains 7™ \ the formula for g(z) can be rewritten as 



/(C) 



(C-^)i?(C)- 



(2-13) 



where 7™ C 5 is a negatively oriented loop around 7™ (which is "pinched" to 7™ in where / is not 
analytic) that does not contain z. Introducing function h = 2g — f we obtain 



h{z) = 



2TTi 



/(C) 
(C-z)i?(C)- 



where z is inside the loop 7m. The endpoint condition (2-111 can now be written as 



or, equivalently. 



h{z) — 0{z ~ a)^ as z ^ a, 



(2-14) 



(2-15) 



(2-16) 



(C-a)i?(C)- 

The latter equation is known as a modulation equation. The function h plays a prominent role in this 



paper. Using the fact that the Cauchy operator for the RHP (2-101 commutes with differentiation, we 
have 



h\z) 



2TTi 



no 



(C-z)i?(C)- 



(2-17) 



where z is inside the loop 7^. 

In order to reduce the RHP (2-5l-(2-6l to the RHP with piece-wise jump matrices, called the model 
RHP, the signs of 5ft.(z) in the upper half-plane should satisfy the following conditions: 
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is negative on both sides of the contour 7„ 



• there exists a continuous contour 7c (complementary arc) in C+ that connects a and so that 
Qh(z) is positive along 7c. Since '^h{z) > on the interval (— oo,/i_), the point /i_ in 7c can be 
replaced by any other point of this interval, or by — oo. 



Note that the first sign requirement, together with (2-101, imply that Qh{z) — along j„i- Since the 
signs of play an important role in the following discussion, we call by "sea" and "land" the regions 

in C+, where is negative and positive respectively. In this language, the complementary arc 7c 

goes on "land" , whereas the main arc 7™ is a "bridge" or a "dam" , surrounded by the sea, see Fig. [3] 

2.2 Reduction to the model RHP 



We start the transformation of the RHP (2-5l-(2-6l by deforming (preserving the orientation) the interval 



(—00, /i_|_) , which is a part of its jump contour, into some contour 7+ in the upper half-plane C+, such 



that e 7+ . Let 7_ be the Schwarz symmetrical image of 7+ . Using the factorization ( 2-8 1 the RHP 



(2-5l-(2-6l can be reduced to and equivalent one where: 



• the right factor of (2-8 1 is the jump matrix on 7 



+ 1 



• the left factor of (2-8 1 is the jump matrix on 7_; 



• the jump matrix on the remaining part of M is unchanged. 

It will be convenient for us to change the orientation of 7+, which causes the change of sign in the off- 
diagonal entry of the corresponding jump matrix. On the interval (/i-i_,oo) we have 5/o(z) < and it 
appears that the jump is exponentially close to the identity jump and hence it is possible to prove that 



it has no bearing on the leading order term of the solution (2-7 1 (as e — > 0: see |7] for the case when /g 



is given by (2-4 1 and [13] for the general case). Therefore the leading order contribution in (2-7 1 comes 
from the contour 7 = 7+07-. In the genus zero case, the contour 7 contains points a, a, which divide it 
into the main arc 7^ (contained between a and a, and the complementary arc 7c = 7\7m. According to 
the sign requirements (2.1 1, the contour 7,„ is uniquely determined as an arc of the level curve 3/i(z) = 



(bridge) that connects /i+ and a, whereas 7c can be deformed arbitrarily "on the land". Because of the 



Schwarz symmetry 2.1 it is sufficient to consider 7 only in the upper half-plane, i.e., it is sufficient to 
consider 7+. 

Having found the branch-point a, the (7- function g(z) and the contour 7,„, we introduce additional 
contours customarily called "lenses" that join a to fi+ on both sides of 7™ (and symmetrically down 
under). These lenses are to be chosen rather freely with the only condition that must be negative 
along them (positive in C_). This condition is guaranteed by (2.1). 
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The two spindle-shaped regions between j„i and the lenses are usually called upper/lower lips (relative 
to the orientation of 7„j. At this point one introduces the auxiliary matrix- valued function Y(z) as follows 



e = 



1 -e-¥''(^) 

1 

1 e-T'H^) 

1 



outside the lips, 
in the upper lip in i 

in the lower lip in C 



The definition of Y(z) in C_ is done respecting the symmetry in Prop. 2.1 namely 
The jumps for the matrix Y{z) are reported in Fig. [S] 



1 -e 
1 



1 
-1 



1 



-e f 



h 1 



1 



(2-18) 



(2-19) 



V / 

Figure 3: The jumps for the RHP for Y. The shaded region is where 3/i < (the "sea"). The blue 
contour is the main arc 7^, the black contour in C+ is the complementary arc 7c and the red contours 
are the lenses. 



The model RHP. In the limit e — > 0, according to the signs (2.1 ), the jump matrices on the comple- 
mentary arc 7c and on the lenses are approaching the identity matrix 1 exponentially fast. Removing 
these contour; 
jump matrix 



these contours from the RHP for Y{z), we will have only one remaining contour 7„i with the constant 
^ ^ on it. This is the model RHP. Calculating the (1,2) entry of the residue at 



-1 
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infinity (see (2-7)) of the solution to the model RHP, one obtains the leading order term of the genus 
zero solution as follows ([7 ) 

qo{x, t, e) = Qaix, t)e-"-^^o ^<^o{s,t)ds . (2-20) 

To justify removing contours with exponentially small jump matrices, one has to calculate the error 
estimates coming from neighborhoods of points a, a and /i+ (for a and these neighborhoods are 
shown as green circles in Fig. |3]). This is often done through the local parametrices. We shall consider the 
construction of the parametrices near these points as already done and known to the reader. We refer to 
the literature |7] |T3] . The only information that we need is that these parametrices allow to approximate 
the exact solution to within an error term £{z) = 1-1- 0{e). 

3 The analysis near the right breaking curve for NLS 

When x, i is on the breaking curve the asymptotic solution to NLS is about to change its genus. Although 
the phenomenon described in this paper is typical for any increase of genus, we will focus on the first 
break, i.e., on the transition from genus from to genus 2. This happens by either: 

• creating a new band (right breaking curve in Fig. |4]); 

• splitting an existing band (left breaking curve in Fig. |4]). 

Let (xo, io) be on the (smooth part of the) breaking curve; this implies that there is 77 G C+ 
3/1(77, xo, to) = , h'{r]{xo,to),xo,to) ^ , C ^ C{xo,tQ) := -ih"{r]{xo,to),xo,to) ^ 0, (3-1) 
where prime denotes derivative with respect to the spectral parameter z; moreover, '^h(z, xo,to) has the 



required signs (2.11 on both sides of main (bands) and complementary (gaps) arcs. We keep notation 



h{z;x,t) for extension of /i(z;xo,to) in a small neighborhood of {xo,to) even if the sign-distribution 



requirements (2.1) are not satisfied. 
Definition 3.1 We denote the critical value 

S{x,t) := h{ri{x,t),x,t), (3-2) 

which satisfies ^S{xo,to) = 0. 

We will now show that S{x,t) is locally a diffeomorphism from the {x,t) plane near the smooth part 
of the breaking curve, into the complex S-plane. From [J, formula (4.46), we have 

h^{z) = R{z) , ht{z) = 2{z + a)R{z) , R{z) = ^/{z ~ a){z - a), (3-3) 
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Figure 4: The typical zero-dispersion phase diagram for a one-hump initial data. Representative at 
different points in the (x, i)-plane are the level-curves of The genus-zero region in the (x, t)-plane 

is the lighter area, while the wedge-shaped area is the genus 2 region. Only the upper-half spectral plane 
is depicted. The shape of the level curves was obtained numerically on a simple example. 

where the square root is chosen to have the behavior R{z) ^ z near infinity (this is the opposite branch 
compared to the one used in 7 ) and analytic off the main arc. Since (by definition of 77) h'(ri) = we 
have that the partial and total derivatives of S{x,t) = h{ri(x,t), x,t) coincide and hence 

= Rir]) , St = 2(7? + 3?(a))i?(7?) . (3-4) 

We also have 

Lemma 3.1 (Lemma 4.5 in [7]) For t > and for any x eM. 

• the function ^hx{z) — 5i?(z) is strictly positive in the upper half-plane except for the finite region 
enclosed between the main arc and the vertical segment [a, a]. 

• The inequality '^ht(z) = 25(z -|- a)R(z) holds for \z\ large enough and to the left of the zero level 
curve of^ht{z) that emanates from a. It changes sign when crossing this curve or the main arc. 

In addition to this we will show that 

Lemma 3.2 The function S{x,t) is a local diffeomorphism for a; G M, t G IR+ with values m C ~ 
except at the point of gradient catastrophe. 
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Proof. The Jacobian of the map S{x, t) is 
Sft Sx S^ 



J{x, t) — det 



= ^Sx'^St - 3? 5*35, = 3 [SxSt] = 23(?7 + a)\R{rj)\'' ^ 

= 23(77) |i?(r;)|2. (3-5) 



Thus the Jacobian does not vanish on the upper half-plane and 5 is a diffeomorphism as long as rj ^ a. 
Q.E.D. 

Theorem 3.1 Assume that on the breaking curve t{x) there is a unique gradient- catastrophe point 
(xoj^o); ??(xo,io) — ce{xQ,tQ). Then, outside of any (small) neighborhood of the gradient- catastrophe 
point (xoj^o); the function S{x,t) maps a strip-like neighborhood of the left breaking curve diffeomorphi- 
cally to a strip around the negative real S-axis with the genus 2 phase mapped to 35 > 0, while it maps 
a corresponding neighborhood of the right breaking curve to a strip around the positive real S-axis, with 
the genus 2 phase mapped to 35 < 0. 

Proof. The full derivative of 5 along the breaking curve t{x) is 

— S{x,t) = -^h{r]{x,t),x,t{x)) = h'{r],x,t)r]x H- hx{ri,x,t) + ht{r],x,t)t'{x) = 

= h^{r], x, t) + htirj, X, t)t'{x) ^ 0. (3-6) 

The vectors hx{z), ht{z) in C are not parallel to each other for any z G C+ \ {a}. But S{x, t) is a real- 
valued function along t{x), so 4^'^S{x,t) ^ 0. Since 5(a;o,io) = S{x,t{x)) — > -l-oo as x — > -l-oo (see 



(3-18), the right branch of the breaking curve is mapped onto the positive real axis of 5. Similarly, the 
left branch of the breaking curve is mapped onto the positive real axis of 5. The required sign distribution 
of 35(2:, i) in the genus two region follows from diagrams at Figure |4] Q.E.D. 



According to Thm. 3.1 the real 5-axis (up to a small neighborhood of the origin) is diffeomorphic to 
the two branches of the first breaking curve. 

We will consider mostly the right-breaking curve, leaving the (not too different) analysis of the left- 
breaking curve to later: the genus 2 region corresponds to negative values of 35. We want to "explore" 
the genus- two region within a scale £|lne| and hence we introduce below a conveniently normalized 
exploration parameter k. 

Definition 3.2 (Exploration parameter) We define the two functions 'd{x,t), K{x,t) by the relation 

2S{x,t) =-d{x,t) +i>ieliie (3-7) 
This equation defines implicitly a, foliation of the a neighborhood of (I'g, t^). The breaking curve is defined 



by the equation x — and "d serves as a parametrization of it, thanks to Lemma 3.2 The scaling in ( 3-7 1 
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means that k parametrizes an infinitesimal strip-like neighborhood of the breaking curve with a width 
of order £ In £. 
We will write 

2ih{z) = -C{z-^f(\ + 0{z-J])) + 2iS, ri^'n{x,t), S -.^ S{x,t), C -.^ C{x,t). (3-8) 
The complementary arc 7c that goes through 77 can be deformed so that - in a finite neighborhood of 77 - 

^{h{z) -S) = Q when z e jc- (3-9) 



Definition 3.3 (Scaling local coordinate) Define 



2i 



(Hz) - s) 



(3-10) 



where the root is chosen so that the complementary arc 5R(ft, — S) — is mapped to the real C^-axis with 
the natural orientation. 



The jump on the gap (complementary arc) that goes through 779 is 



1 

_e¥M^) 1 



= 11 



1 

— e + e 1 



2 2i!R(S) 



,2 2i 

-e-^e-'^ — 



=: y_ 



1 



-e "e 



(3-11) 

Note that the oscillatory term e*''/^ is not problematic because it has constant modulus, but the term 
diverges for x > Q. 

3.1 Long time asymptotics of S{ri) along the breaking curve 

Let us recall that the scattering data /o(z) is analytic in C-|_ (solitonless case) and that the corresponding 
w{z) = Q/o(z)sign(z — /i+), restricted to z e M, is admissible in the sense of [13]. We recall that the gist of 
the notion of admissibility for w is that |»'o(2;)| has exponential growth on the interval (fi^,fi^) as e — s- 0, 
exponential decay outside and w is piece-wise continuously differentiable with w'(/i±) 7^ 0. As 

it was mentioned above, the condition of the absence of solitons is not essential if the contour 7 of the 
matrix RHP can be chosen so that it does not intersect singularities of f(){z). 



Under the above assumptions and using (2-17), (2-141, the system of equations for the point rj = ri{x, t) 
and the breaking curve t = tQ^x) can be written as 



2t(77 + a) + x + 



wiOdC 



1 



ic-v)m)\ 



= 



= 



(3-12) 
(3-13) 
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According to [13], Sect. 3.4, this system implies that the right branch of the breaking curve has the 
asymptotics 



t{x) 



1 



X + 0(1) , X +00 



and correspondingly 



^ =(/._ + or l))+^(^^[l + o(l)] 



t oo. 



(3-14) 



(3-15) 



Moreover it was shown in Lemma 3.6 [l3j that the branch-point a{x, t) approaches /i+ as a;, t oo along 
the right branch of the breaking curve according to the following asymptotic formula 

2,rt(M+-M_) 



~ja = e 



-+o(i) 



The two later equations, together with Proposition 4.3 from [TS], applied to 



Hv) = Riv) 



2t{r] -I- a) -I- X + - 



^J-oo {C-il)\R{0\ 



yield 



S{'q) - At^i^in^ - fi+) 
as ?7 oo along the right branch of the breaking curve. Using 



h"iv) = 



2m A (c - vrmcY 



/'(C) 



-dc 



(3-16) 

(3-17) 
(3-18) 

(3-19) 
(3-20) 



we can, in a similar way, calculate 

C{x,t) 4it 

along the right branch of the breaking curve. Similar asymptotic analysis can be done for the left branch 
of the breaking curve. 



In the case of fo{z) given by (2-4l, using the explicit expression for h"{z) given in [Tj, Section 4.5, we 
obtain 

. n — n IT 

(3-21) 



S{t)) = ~2m(r]) h + a] + sinh"! ^ - z| 



4 First approximation 

4.1 Outer parametrix 

The construction of the outer parametrix, namely, the solution of the model problem in Section |2.2| is 
accomplished along the lines of f5^. The outer parametrix (aka model problem) is the unique solution of 
the RHP 



*+(z) = v|/_(z) 



1 

-1 



on (a, a) 



^{z) = 0{{z-a)-^) ^{z) = 0{{z -a)--^) 
«'(z) = 1 + ©(z^i), z^oo. 



(4-1) 

(4-2) 
(4-3) 
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^'(z) is explicitly given by (4-13 1 below. Following the ideas in [9] we introduce the so-called Schlesinger 
transformation of the model problem 

Problem 4.1 (Schlesinger-chain of model problems) Find a matrix ^'^(z) analytic off [a, a] U 
{77,77}, rj ^ [a,cF], such that 



1 

-1 



■^^{z) = 0{{z^a)-i) 
vl'^(z) = l + 0(z-i) , 
*,(z) = 0(l)(z-77)-^--^ , 
M/,(z) = 0(l)(z-77)^-^ , 



077 {a, a) 

Vl/,,(z)=0((z-«)-3) 



> 00 

z rj 
z fj. 



(4-4) 

(4-5) 
(4-6) 
(4-7) 
(4-8) 



We have formulated the problem assuming that there is only one main arc (cut) but the modifications in 
the case of multiple arcs is straightforward (see [9J); the main point is that we have modified the model 



problem only by adding the behavior (4-7 1 and (4-8 1 near the points where a new cut is about to emerge. 



A- 



a, 



(4-9) 



Uniformization of the plane without the cut. Consider the complex z-plane with a slit on the 
segment [a, a]. It is elementary to show that the map 

a — a/, 1\ a + a b 

where a = a+ib, maps univalently and biholomorphically the 
outside of the unit disk in the A-plane {|A| > 1} to C \ [a, a], 
with the unit A-circle being mapped onto the segment two- 
to-one. The points A = ±7 are mapped to a, a respectively. 
If we choose another (simple) curve joining a to a then its 
counter-image in the A-plane will be some deformed curve 
passing through A = ±7 and encircling the A-origin. Note 
the symmetry 

z(-A-i) = z(A). (4-10) 

Therefore the outside and inside of |A| = 1 are conformally 
mapped to two copies of z £ C \ [a,a]. We will refer to 
the region |A| > 1 as the physical sheet, as opposed to the 
unphysical sheet |A| < 1. 

The points z = r],rj have unique pre-images Ag, Aq on the 
physical sheet 

^(^o) = V , z(Xo) = T] , 

We then have 



■ \) 


A 

1 i 








/ * ' 

/ 1 ' 














A-planc 


\ 1 1 






• \] 




i 








a 














z;-plane 


a 





Figure 5: The uniformizing plane: the 

shaded region is the unphysical sheet. 
|Ao| > 1. (4-11) 
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Proposition 4.1 The solution of the RHP (4--1) exists (unique) for any K ^'L and it is given by 



Kc73 



AAo + 1 



A_-Ao 
AAo + 1 



Ka3 



\=X{z) 



(4-12) 



where A(z) is the determination on the physical sheet o/ z G C \ [a, a] and Ag (defined in (4^-11)) is the 
counter-image of rj on the physical sheet. The matrix ^'^ is given by 



1 


—i 


-1 ■ 




i 


1 




/ \ 0'2/4 


2 


1 


i 




-1 


— i 







Proof. It is well known that the matrix "ifg solves the RHP 

for K — 0. Since the interchange of sheets corresponds to the map A 1 

the Szego function 

D{z) 



AAo + 1 A - Ao 



A - Ao AAo + 1 



(4-13) 

' ^ one immediately verifies that 

(4-14) 



satisfies D+(z)£'_(z) = 1 on the cut, and limz_>oo D(z) = =A. This is the reason for the left normalization 
in eq. [4121 Q-E.D. 

In terms of A = A(z), the matrix is given by 

1 



(Al - iCTa) = 



A -1 
1 A 



(4-15) 



Remark 4.1 Since = (A — the uniformizing variable A is related with the hyperbolic variable v, 
defined by sinhu = from Sect. 4-5, through \ — e" . Thus, 

z — a + R{z) 



X{z) = e" = 

where the square root R is chosen so that A(z) is on the physical sheet. 
Lemma 4.1 The matrices ^'^ are related by 



(4-16) 



^,^,{z) = RJz)^fJz) , RJz) = l- 



C 



K 



C 



K 



Qk _ Qk 
z — rj z — rj 



(4-17) 



z — rj z — rj 

Proof. The matrices k have the same jumps and the same asymptotic behavior at the branch- 
points a, a and at infinity. Thus Rk ■= ^'^^-^^^'^^ is a holomorphic matrix (without jumps) with possibly 
isolated singularities only at ry,^. The asymptotic behaviours of and at these two points imply 

that Rk may have at most simple poles there. The same reasoning can be applied to ^'^^k+ij whence 
the formula for i?^^. Q.E.D. 

It is possible to give explicit formulae for the matrices Ck^Ck^QktQk in terms of the columns of 
^ K but they are of no immediate interest. 
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4.1.1 Reality 

The matrices ^'^(z) can be written as 



Ak(z),Bk(z) 



[z-ri) 



(4-18) 

where the columns Ak,^k are analytic a,t z = rj and, respectively, the columns AK,JiK are analytic at 

z = f]. 

Note that, since det^/f = 1, we have 



det 



Ak{z),Bk{z) 



= 1 = det 



Ak{z),Bk{z) 



Due to the symmetry of Prop. |2.1| we have 

'i>K{z){^K{z))^ = 1 , 

which implies ^K{z){^K{z)y = 1 or 

[Ak{z),Bk{z)]= (JAk{z),Bk{z) 

This is promptly shown to translate into the relations 

1 



^k{z) = ( ^ Q 1 Bk{z) , Bk{z) = ^ 



-1 



Ak{z) 



(4-19) 



(4-20) 



(4-21) 



(4-22) 



$k(z) = CT2«'k(z)ct2, 
which is equivalent to ^K{z){^K{z)y = 1- Similar considerations lead to 



-i 

1 



$'^(z) = cr2<^>'K{z)a2 , (72 ■■ 

We will denote by Aj^ the value Axirj), and similarly 

Ak := AKiv) , 6^(7?) , A'j, := Ai^(z)U^^ , B'^ := B'^(z)|,^^ , 

:= AKiv) , B^f Bk(??) , A'x := A'K(^)U=rj , B'^ B';^(z)U=^ . 

This implies the following relations, which will be of use at a later stage in the paper 



det [A'^, Aa-] = det 



det 



, Bk 



B^^Ba > 



det [Ba,Ba] det 



A'a, Ak 



det 



Aa, B^- 



AT^Aa > 



These quantities are computed explicitly in App. \K\ 



(4-23) 



(4-24) 



(4-25) 
(4-26) 



(4-27) 
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4.2 Local parametrices near ri,ri and first approximation to the solution 



We only consider the point 77, with similar considerations being repeated for rj. The local parametrix 
must solve the following mixed RHP where the adjective "mixed" means that the requirements are for 
both columns and rows: 



H, + (z)-H,._(z) 
H^(z) = 1 + 0(1) , zee 



(4-28) 

(4-29) 
(4-30) 



where o(l) above means an infinitesimal in e, uniform w.r.t. z when z belongs to the boundary of the 
disk. When necessary to specify the values of the parameters >c, z9 we will write 

H,+{z)^H,+{z;>c,d) . (4-31) 
4.3 Approximate solution and error analysis 

Suppose that we are able to solve the mixed RHP in ( |4-30[ ) and define the following matrix 

{^'/f (z) outside of the local disks , 
^'a'(2)Hx(z) within the disks at 77, , (4-32) 
'i' k{z)A{z) within the disks at Q;,a , 
where A{z) is the usual Airy-parametrix [H]. Note that the matrix ^'x(z)Hif(z) is bounded in a 



matrix Y(z) we take their ratio 



neighborhood of 77 because the singular behavior of the columns of ( |4-7| is exactly canceled by the 
singular behavior (4-291 of the rows of H^. In order to evaluate how close the matrix Y{z) is to the 

£{z) ■.^Y{z)Y{z)-^ 



and consider the RHP it solves: 

£{z) ~ l-l-O(z-i), 



£+{z) = £^{z)W{z) 



(4-33) 



(4-34) 



where M^(z) is a matrix supported on the contours indicated in Fig. [6] In 
particular, on the boundary of the disk D centered at 77 we have 




£+{z)^£_{z)^k^k{z)^k\z) ^ ze9D. (4-35) Figure 6: The con- 

tours for the RHP 

Since are uniformly bounded for z G 9D as e — > (in fact they are independent of the error term (in 
of e), it follows that the jump of £ on 9D is 1 + o(l), where this is the same '^+)- 
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infinitesimal that appears in eq. (4-30) (and we will promptly see below that this desired infinitesimal 
is not such when x is a, half-integer, thus requiring a more refined approach). The usual theorem about 
Riemann-Hilbert problems with small-norm jumps imply that £{z) is uniformly close to the identity to 
within the same infinitesimal. 



We will see in Section 4.4 that it is possible to solve the local problem for H^f formulated above with 
o(l) = C'(£5-l"--f^l) and hence: 

• the integer K must be chosen as the nearest positive integer to x\ 

• it is impossible to solve the problem as formulated for a half-integer >r = ^ + A/", where A/" G N. 

We will show below how to solve the problem as formulated, why it is impossible to solve in those 
exceptional situation and how to circumvent the problem. 



4.4 Solution of the local RHP for Hk (4-28). 



Let Hi{s) be the monic Hermite polynomials of degree i: 
We define the local parametrix as 



ds = -^^Ow ^ifdw , 



H.is) (-2)"V — e- 
ds* 



(4-36) 



Hk(C) 



-1 



K-l 



{s)e- 



-2i7r 



-c 



K 



Jr s - C 
1 r HK{s)e-'^ds 

2i7r s - C 



Hk -liC) 



HKiO 



2iTT 



HK{s)e-"' ds 



H 
Hk{0 



K-l 



(0 



5-C 

Boundedness. The matrix ^'^Ha' is bounded in 



c 



K 



(4-37) 



(4-38) 



because of the (^^'^'^ factor in (4-37) cancels out 



the singularity of k (4-7 1. 



the zooming coordinate C, scales like Therefore the 



0{e 



), ze 



(4-39) 



Error estimate. Recall that for z £ 
matrix on the boundary of the disk D has the following behavior uniformly w.r.t. z 

_ r l + 0(e) ©(£3+--^ 
" [ 0(e5-^+^^') l + 0(e) 

Note that the diagonal estimates are 1 -I- ©(e) and not just 1 -I- 0{y/e) due to the fact that the Hermite 
polynomials are even/odd functions according to the parity of their degrees. For example, for the (1,2) 
element in (|4-37[): Hk-i{0 = C^^^i'^ + C(C"^)) and on dB we have C = 0{l/^/e) and 



2i7rC^e~'-£^ 



2mC 



2K-1, 



-AH 



■ 



(4-40) 
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Similarly, for the (2, 1) element, using the orthogonality, we have that the Cauchy transform is ©(C ^ ^) 
and hence 



C-'^e*f£- 
2z7r 



s — 4 zin 



(4-41) 



We immediately see that -in order to meet the requirement about the decay on the boundary of the disk 



(4-30)- we must have 



Proposition 4.2 The integer K must he the closest integer to x. 



The important observation is that if x e I + ^ then the error term in ( 4-39 1 does not tend to zero (it 



is 0(1)). It is understandable as these values separate regimes where the value of K jumps by one unit 
and the whole strong asymptotic must changes its form. 



5 Improved approximation 

5.1 Next order improved local parametrix 

Inspecting directly and a bit closer the expression ( 4-37[ ) we see that 

-2i7rC 



H 



K 



1 + 0{e) 



0{e 



+ 0(ef— +^) 



1 + 0{e) 



(5-1) 

and the terms that are responsible for the loss in the error term (4-39) are the off-diagonal ones. In 
particular: 

• for K < K < K + ^2 the leading term in the error is the lower triangular one; 

• for K — ^ < K < K the leading term in the error is the upper triangular one. 
For this reason we re-define 

Definition 5.1 The local parametrix H^, >c G is defined as 



1 



2iTlQ 



= < 



H 



K 



K <x<K 



K 



(5-2) 



2i7rC 



fJ'K-1 
1 



H 



K 
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Parametrix on the disk around 77. Because of the symmetry of Prop. |2.1| we have 

In the formulae below, we will use C(z) :— C{z). 
5.2 Improved outer parametrix 



(5-3) 



Corresponding to the definition of H^^^ (Def. 5.11 we need to define the outer parametrix with the 
requirements 



1 

-1 



^';^(z)+ = ^'^(z)_ 

*^(z) = 0{{z - a)-i) ^>^{z) = 0{{z - a)"4), 
*^(z) = l + 0(z-i), 
*^(z)H^(C) =0(1) , zeD, 
■^^{z)n^{C) = 0{l) , z e D. 

In correspondence with Def. |5.1| we will pose the following Ansatz for ^f^^. 



(5-4) 

(5-5) 
(5-6) 
(5-7) 
(5-8) 



Theorem 5.1 The solution of the Riemann-Hilbert problem (5-4-)-(5-8) always exists and is given by 



where 

GkA*) = 



TkA^)^k{z) : 
UkAz)'^k{z) 

uk det[»,B;<-] 



1 - 




Gk,>c 




z — 77 


Z-T] 




(^) 












Z — Tj 


z — fj 



^k{z) K <k<K+\ 

>f = A' e N 

-\^k{z) K~\<x<K, 



(5-9) 



l + UKdet[B^,BK] 



l + ?lKdet[B'^,BK] Bk- 



uk||B"||2 



4a()))2 



23(r/) 



(72B^S-,10) 



det[Aif,»] 



l-^ifdet[A'^,Ax] 



49'(?))2 



— (^(1 -Z;^det[Ai,, A,,']) Ak ~ ^^^^|^^2Ak^,11) 



2^+ii7r 

The proof can be found in App. iB 



2i7re-''?e^--^+5 



(5-12) 
(5-13) 
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Remark 5.1 The idea behind this Ansatz is the following: the outer parametrices k o,nd ^'a'+i are 



related by a finite Schlesinger transformation of the form (4--11) and we know that they provide the correct 



-\ respectively. The discontinuity in this 



leading behavior for x E {K — ^, K + ^) and x E {K + ^, K , ^ 
leading-order asymptotic behavior occurs at x = K + ^, where the error term as per (4-39) becomes of 



order 1; the idea is thus to introduce some sort of interpolation between the two outer parametrices of the 
same form of a finite Schlesinger transformation, letting the coefficients matrices be dependent on k in 
such a way that their order becomes I as x crosses the critical value. 



Suppose K > K + ^ in the formula (5-10) for Gk.>c {Gk.>c)', then one sees that uk diverges as e — s- 0. 

{GK.>t — > Ck), while ior x < K 



K 



The matrix Gk,3< {Gk,x)i however, has a regular limit Gk,. 
we have Fk,>c — > Qk , {Fk,>c ~> Qk) where Ck, Ck, Qk, Qk have been introduced in (4-17); indeed the 



defining equations for G,Gk,>c turn into those of Ck,Gk (ditto for Fk.:k, Fk.}< versus Qk,Qk)- This 
means that Gk.>c,Fk.>c are some sort of interpolation between Schlesinger transformations. 



Continuity at x = K + ^. We must check that the matrix ^ as defined in eq. (5-9) is continuous 



as X = K + at least within the error ©(-y/e). In fact we are going to see below (Prop. 5.1) that this 
continuity holds exactly, namely. 



lim = lim ^ 



or, equivalently, 

/ 1 Gk.. 



G 



K.: 



z — rj z — ?; 



Fk, 



Fk+1 



Z — T] 



z — rj 



"^K+iiz). 



(5-14) 



(5-15) 



Even more is true as eq. (5-15) holds for any value of h. To show this we recall that ^k+i = FLk^ k 



(as per eq. 4-17). Thus, we prove the following 



Proposition 5.1 For any value of >f G M, the matrices Gk,3<,Gk,k, Fk+i.^, Fk+i,>c are related with 
the Gk,Ck+i by 



G 



K.: 



G 



K.: 



Fk-\ 



Fk 



Gk Ck 



z ~ rj z — rj 
In particular, we have 



z — rj 



z — rj 



z — Tj z — rj 



F, 



Fu 



Ck + C}s 



(5-16) 



(5-17) 



The proof can be found in Appendix |B] 

Thanks to Proposition |5.1| we conclude also that in Theorem |5.1 



Tk.^{z)'^ k{z) = Uk+i,>,{z)'^k+i{z) , "ix, z, s, 
which allows us to dispose of the necessity for the Fk altogether and rewrite Theorem |5.1| as 



(5-18) 



26 



Corollary 5.1 The solution to the RHP (5-4)-(5-8) given by Theorem 5.1 can he equivalently written as 
^^(z) =[l- ^ - ^ ) ^k{z) , K = max(0, W ) , (5-19) 



z — rj z — Tj 



with the matrices being defined in (5-10) 



We remind the reader that k (with integer subscript) had been defined independently in Proposition 

Remark 5.2 We also remark that the new definition in Cor. \5.1\ is not continuous in x at the integers 
(but continuous at the half-integer). However the discontinuity is to within 0{^/e). More precisely, 



lim *^(z) = ( 1 + 0(V^) ( ^— + 
<^K- \ \z~ri z — T] 



lim 'i'Jz) 



(5-20) 



with the 0{y/e) factor being independent of z. This discontinuity is of no concern because the correction 
to the final approximation is within our overall error estimate. 



6 The left breaking curve 

The situation at the left breaking curve (Fig. |4]) is at first sight quite different because it is the main arc 
that it is splitting into two, rather than the complementary arc. However, as we will presently see, the 
phenomenon is essentially identical: the informal reason is that the choice of main and complementary 
arcs in the construction of the function h is somewhat conventional, depending on the overall sign of h at 
infinity and the placement of the cuts of the ^-function. The difference is nevertheless sufficiently stark 
to grant a separate analysis. 

At the exact point of break the function h near 
the point rj (where the arc is about to split) can 
be written 



2ih{z) 



^C{z-r,Y{l + 0{z-ri)), 
C{z-iqf{l + 0{z~ii)), 



z e Oi, 
z e Do, 
(6-1) 

where the region Di is the region (within the lo- 
cal disk) to the right of the two main arcs and Dq 
its complement. If we let evolve the level curves 
according to the modulation equations while re- 
specting the branch-cuts we would have the situa- 
tion depicted on Fig. [8] with two pairs of close cuts 
running off to infinity (i.e. a new main arc appears 
to compensate the split of the original one). 




Figure 8: The main arc has split into two and we are 
forced to introduce a new main arc on the zero-level 
set of 3(/i). 



27 



Figure 7: The typical shape of the main, complementary arcs and lens contours at the precise point on 
the left breaking curve. 



We then must re-define the function h by introducing a small main arc within the local disk D 
according to Fig. [9] 

The difference between the two pictures on Fig. [8] and Fig. [9] is that we have changed the sign of h 
within the region that has changed color; the result is a function that we still denote by h(z), which is 
harmonic analytic on the upper complex plane minus the main arc (marked by the thick line). In the 
sequel we will denote by h{z) the function defined this way: it has the property that it is now analytic 
in a small neighborhood of z = rj and behaves 

2ih{z) = -C{z - rif{l + 0{z - rf)) + 2iS, (6-2) 
C^C{x,t), ?7 = 77(a;,i), S -.^ S{x,t) . (6-3) 

The arc made by the two lenses can be deformed to pass through the level-curve 3?(/i — S*) = and the 
zooming coordinate C, can be defined as the conformal change of coordinate 

r f(M^)-5), zGDi, 
C = < (6-4) 
i -^-1 {h{z) - S) , zGDo. 

We note that now we have 5(5') > and hence it is convenient to parametrize it by 

S -.^id + iixelne, (6-5) 
where -contrary to the case of the right breaking curve- negative values of k correspond to the explo- 



ration into the genus-two region. In Fig. 10 we have shown the jumps and contours for the solution Y{z) 

1 
-1 



in the local disk D. The model problem solves the usual jump = ( ^ ^ 1, 1 
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Figure 9: [On the left] The main arc has spht into two and we are modifying the function h by introducing 
a small "artificial" main arc connecting the two pieces of the broken one. 

Figure 10: [On the right] The exact jumps of Y in terms of the conformal coordinate C near the point rj. 

6.1 Outer parametrix 

The situation from this point on is quite similar to the previous; the outer parametrix ^'^ will be taken 
as defined by the same Riemann-Hilbert Problem |4.1[ with the only difference that the jump condition 
now runs on the piecewise-smooth arc. Note that the point r/ (by our construction) lies to the right of the 
main arc; worth mentioning is also that this time we will have to use the formula for \E'^ for the closest 
negative integer to >c. The formulas for "ifj are identical. 

6.2 Local parametrix 

Due to the different shape of the jumps within the local disk D, the relationship between the present case 
and the previous ones are less immediate than for the model problem. 

Problem 6.1 The local parametrix is a matrix-valued function Vj^C.) — V x^d) of the zooming 
coordinate C, satisfying the following mixed RHP G — Nj. 



1 






(6-6) 







) 1 

1 



Q G main arc. 



(6-7) 



(6-8) 
(6-9) 
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We anticipate that we will need this problem for negative integer K. If we disregard the condition (6-7 1 



then this problem is essentially (up to the triangularity of the jump) the same problem as for H^, ; more 
precisely, let K — N and x e — M+ and consider 



(6-10) 



is the solution of (4-281, (4-29), (4-301 with K = - J and ^tr, replaced by -m:,--^. Then the solution of 
the RHP [O is 



1 

1 



H, ,,(C;->^,-^) 



1 

1 



VAC,^,^) :=e^^^H,_,,(C;-x,-^)e-^' 



(6-11) 
(6-12) 



It is promptly seen that V , solves the jump on the real axis and -on the boundary between Do,ID)i- is 



satisfies the condition (6-7) 



It should be apparent now that the sequel of the analysis is quite parallel to the previous case; the 
parametrices ^ and will be constructed by the same methods. The explicit formulas for the matrices 
Gj,Gj etc. have only minor cosmetic differences. 

7 Corrections to the NLS solution 



Theorem 7.1 The solution of the NLS (1-1) defined by the scattering data fi){z) behaves as 



as e — > 0. Here {x,t) are such that 



S(S'(x,t)) = -elne , >ceR, 



(7-2) 



where K = max(0, \_>c\), e^"^ :— (for Ao see (4--11)), Gk,>!,Gk,>c have been defined in {5-10). The 



term {Gk. 

2(G'_R-,jf -I- GkA^'^ 



Gk,>c)i2 has the following behavior 

( 0(£3+^"-"), 



-b{x,t) 



(7-3) 



In particular the correction is of the order 0(1) one as e — > only for k E ^ -\- N. 
The proof can be found in Appendix [B| 

Remark 7.1 Considering x as a parameter in a direction transversal to the breaking curve, we notice 
that centers of adjacent ripples are at unit x-distance from each other. On the other hand, according to 



{7-3), the "support" of each ripple is of the order 0{^^) in terms of the exploration parameter x. Since 



the order of x is ©(elne), this means that the support of each ripple is of order 0(e), exactly as in the 
bulk of the genus 2 region. However their spacing is of the same order as x, so that they look like ranges 
of mountains/valleys separated by vast flatlands (see Fig. 
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While Theorem |7.1| contains the substantial quahtative information about the relevant features of the 
behavior, the reader may wish for a slightly more concrete formula. It is possible to write the correction 
term in full detail and -for example- use some computer algebra software to achieve the plot in Fig. [2j 
We restate the theorem with some more explicit expressions. 



Theorem 7.2 The solution of the NLS (1-1) defined by the scattering data fo{z) behaves as 



qix, t) = -bix, t)e^+'^'<^ [l - ^^^^ ^^^^ I +u,B,/ + B,,2 j + 0{^e)^ , (7-4) 



wh' 



ere 



- — ; , Bk,2- ,o : (7-5) 

1^1' z63(Ao)(l-HAg)2' ^"'> 
= -I ^' ^, , exp ( £3--+^ , (7-7) 

'''' ^ll + A^lHl + Ao 1 + ^0 j- ^'"'^ 

The proof can be found in App. |B] 

7.1 Discussion 

Universality. The correction depends only on "^k, V ^'^d the constant C that appears in uk', in 
particular, k depends only on the branch-point a and the integer K which counts the number of first 
ripples. Thus, universality here is taken to mean that no details of fo{z), that defines a particular solution 



to the NLS ( 1-1 1 , enter into the formula. 

Since the gradient of ^^(x, t) does not vanish, the natural scale of modulation of q(a;, t) in the direction 
transversal to the breaking curve is eln iUVSiSIl and >c serves as exploration parameter in this scale. 

Consider the case k — K + ^ E ^ + N, i.e. on the "crest" of the i^T-th ripple; for this value of the 
parameter the correction to q{x,t) is of order 0(1). In this case uk is a fast oscillatory function: if 
we fix a small e then the correction term to the amplitude has a slow modulation in (x, t) due to the 
dependence of C, ry, a etc. and a fast change along the direction parallel to the breaking curve due to the 
fast oscillation of u„. 
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Large time limit of the profile. According to [15], [13], together with (4-161, we have the following 
asymptotics along the right breaking curve as t — > oo: 

a = a + ib ^ fi+ + ie , (7-10) 

rJ^^,_+^^, (7-11) 

Ao = A(^) ^ ^^to^ + (7-12) 

^ - -'^^Wrt-v ' "^(^) - » (/o(^)+) ^^9n{x - ) (7-13) 



Then 



Bl,^e'^^^^'^(^^y\ (7-14) 



^-u;'(pi_) 

whereas det[B^,B/f] and Bi^i are decaying exponentially as t oo. Thus 

iK+2 



rK=^K--l+\uK?[^^,] (7-15) 



and, according to (3-20) and (7-7 1, 

UK ^ e 4e'^ ?^-r • (7-16) 

In order to gain some insight on the profile and relative height of the ripples for large time we take >f a half- 
integer (similar consideration hold in a double-scaling approach as long as the expression £^+i~^t^+k 
approaches infinity). Then we obtain {x = K + ^) 

„ oK If \2K+2 

{Gk^k^.-Gk^k^.),, - -V^e-^z-e--'>e'?^^^|^^^ . (7-17) 

In the particular example, considered in the paper for numerical computations, w'{ii-) — ^tt. Notice 
that the ripples has algebraic decay as t — > oo whereas the amplitude of the corresponding genus zero 
solution decays exponentially fast. Thus, F^ippie is the main contribution to the solution along the ripples 
for large t. Expression (7-171 also shows that the rate of decay for large t increases with the number K; 
however, it may be offset by the Kl growth of the coefficient in front of t^'^^ i . That is consistent with 
0{t^^) magnitude of the solution in the bulk of the genus two region, see [TSj. 

Remark 7.2 Similar results are true for transition through any breaking curve separating regions of 
genus 2N and 2N + 2. 

Remark 7.3 This type of scales are to be expected from the similarities with random matrix models J 161/ . 
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A Explicit formulae 



Using (4-131, the expression for <i>_R-(z) can be written as 

/A \'^'^^ 



To express it in terms of A, we calculate 



{z-7^)D{z) = 



&(A-Ao)(l + AAo)- 
2AAo(l + AAo) 



so that (A.ll becomes 



1 /Ao 



A -1 
1 A 



/ b(A-Ao)(l + AAo)- 
V 2AAo(l + AAo) 



The quantities involved in (4-271 require the computation of (f>'^(z), which is given by 



Aq\ I z — a 



2i?2(z) 



Ao / \ z — a 
The expression in the square brackets is 

2A2 



where 
Thus, 



h{l + A2) 
2Ao 



+ KcF:i (ln(z-77)-|-lnD(z))' 



-0-2 + Kv{\)(Tz 



1 + A' 



Ao 



1-FAAo 1-hAAo A-Ao A 



(A.l) 
(A.2) 

(A.3) 

(A.4) 

(A.5) 
(A.6) 



2A2 



A(A-Ao)(l + AAo)- 
V 2AAo(l + AAo) 



(A.7) 



B Proofs 

B.l Proof of Prop. [SJ! 



Equation (5-171 follows immediately from (5-161 by comparing the behavior of both sides of the equation 



as z oo. To verify this we recall that the matrices Gk,x, Gk,x, Fk+i,>c, Fk+i^^k, Ck, Ck are uniquely 



characterized by the equations ( 5-7 5-8 1 that can be written as 



Z — f] Z — T) 



[Ak{z),Bk{z)] 



Ak{z),Bk{z) 



1 

UK 



z — rj^ 

UK 



z ~ rj 
1 



0(1) z^-q 
0(1), z^% 



(B.l) 
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1 _ Ee±^ _ ^|±i^) [Ak+i(z),Bk+i(z)] 



Z — Tj 



z-r-i 



z-rj 



Ak+i(z),Ba'+i(z) 



1 

z — ?7 
1 
1 



1 



= 0(1) z^r/ 
= 0(1), 



z — J] 



Ak{z),Bk{z) = Ak+i{z),Bk+i{z) {z-fjy^. 



We note that from ( 5-13|B.3 1 it follows 

UK^K+l = 1 = Uxf-K+l- 

So the check consists in verifying that the matrix 

\ z — Tj z ^ rj I \ 



Ck Ck 



solves 



M{z) [Ak{z),Bk{z)] 



M{z) Ak{z),Bk{z) 



1 

UK 



Z ^ T] 
^ UK 

z — rj 
1 



— rj z — rj 
= 0{l) 

= o{i). 



Indeed we have 

M{z) [Ak{z),Bk{z)] 



1 

UK 



Z ~ f] 

z — rj z ~fj j 



Ck Ck 



z — rj z — rj 



[AK(z),B/f(; 



1 

UK 



[From ([R4|] = 1 - 



Fk+i.>c Fk+1,} 



[From ( |R2| ] = 0(1) 
= 0(1) 



z — rj z — rj 



z — rj 

[Ak+i{z),Bk+i{z)] {z-Tj)-''' 



1 

UK 



. z — rj 



z — rj 
1 



(z - rjY 



1 

UK ^ 

z — rj 



1 — ^if + lMK 

w_R- (z-rj). 



(B.2) 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



(B.7) 



(B.8) 



This last expression is analytic because Ik+iUk = 1- This proves the identity (5-161. One should repeat 
the verification near r] but it is completely parallel. This proves the continuity. Q.E.D. 
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B.2 Proof of Thm. [O 

The existence of the solution will be implied by the formula and by noticing that all denominators are 



strictly positive. The symmetries (5-12) are obtained from the symmetries (4-231, (4-24) 



The conditions that determine uniquely Gk.xiGk.}<t Fk.:ki Fk.}< are to be read off (5-7 5-8): we will 



show the details only for the pair Gk,>c, Gk,k, the computation for Fk.^c, Fk.:^ being completely similar. 



Case K < X < K + ^. The conditions ( 5-7 5-8 1 translate into 



Gk.>c Gk.> 



z — 1] z ~ rj 
Gk,>c Gk^ 



\ Z — 71 Z — 7] 

This can be rewritten 

/ , Gk.>€ Gk.x 



[AK{z),BKiz)]{z-v)-''^' 
Ak{z),Bk{z)] (z-^)^"^^ 



1 



2i7rC2^+ie-*?£- 

2i7rC2^+ie'fe> 
1 



(2-77)^-^'=0(l), (B.9) 
(z-^)--^^^ = C'(1)(B.10) 



Z — J] Z — T] 

Gk,>c Gk,^ 



z — J] z — rj 



[Kk{z),Bk{z)] 



Ak{z),'Bk{z) 



z — rj 



2K 



-K 



1 



z-rj\ fiK 



1 



(B.ll) 



0(1). (B.12) 



From (3-8) it follows that 

^(^VC{z-rj){l + 0{z-rj)) , Vs( ^ ^^{z ~ f])a + Oiz ~ fj)). 



In particular, we have 



lim 



z ~ rj 1 



lim 



z — rj 1 



and that one can replace the above condition by the following 



1-^-^1 [Ak{z),Bk{z)] 



z — rj z — rj 

Gk,>€ Gk 
z — rj z — fj j 



1 

UK 



Ak{z),Bk{z) 



. z — rj 

^ UK 

z — fj 
1 



uk 



We remark that 



Uk = -Uk ■ 



(B.13) 
(B.14) 



(B.15) 
(B.16) 

(B.17) 
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This implies the following system of equations 



7] — fj 



= 



(B.18) 



-Gk,^'Qk + uk{ Ax — — — Gk,>c^'k 



where -for brevity- Ak ■— Aa'(?7), := A'j^{ri) etc. The solution is 

UKdet[»,BK]{v - 



GkA<) 



GkA') 



(1 + UK det[B'^, B^]) [I - UK det[A^, A^] j (77 - ?/)2 + ukUk det[AA', B^]^ 

( (^ - A i^^i A A T3 , ^^i^-detpK, Ak] -r- 
■k\\\-uk det[AK, Aif] j B^ H {rj-r)) 

-UKdet[»,AK]{ri - vif 



(1 + UK det[B^, Br-]) [l - UK det[A^, A^] j (ry - r/)2 + ukuk det[A/f , Bk] 

1/1 , A ^\T}' Tj RA w/f det[Aif,BK] 
X + det[B'x,Bif])Aif ^ 'Bk 



iv-v) 



where we have used that 



det[AK{z),BK{z)] = det[AK{z),BKiz)] = 1. 



X (B.19) 

, (B.20) 
X (B.21) 

) , (B.22) 
(B.23) 



Using the reality relations (4-271 and (B.17) we can rewrite them in the slightly simpler form which 



is in the statement of the theorem. 

Case K — ^ < k < K. The derivation is completely analog. In particular we will see below that the 
formulae for Fk,k, Fk,3< are not necessary for our purposes if we have those for Gk,^, GK,>e- Q.E.D. 

B.3 Proof of Thm. PtTTI 



According to (2-7 1, the solution to NLS is recovered from the solution T{z) of the onyma/ RHP as 

q{x,t,e) = -2 lim zVi2{z). (B.24) 

z — yoo 

Away from the cuts the approximation is such that 

g-f ff(oo)<r3p(^-)gf = £{z)-^^{z) , (B.25) 
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where £{z) is the error term matrix which is (miiformly) £{z) = 1 + 0{y/e) and, according to (2-121, 

g{z;x,t) ^ g{oo;x,t) + 0{z~'^) , (j)eR. (B.26) 

Thus we have 

q{x,t,e) = -2e^a{oo;x,t) ^.^ z^^i2{z). (B.27) 

z — >oo ' 

According to Thm. |5.1[ Cor. |5.1|and Proposition |4.1| 



Aq 
A(z) 



K = max(0, 



(B.28) 



Then 



■2 lim zi^-^lis = -2 



2K 



lim z[*o]i2 + 2(Gi<-,^ + G/f,;^)i2 



&(a;,<)e^^^'^+2(GK,.-GK,,.)i2 



(B.29) 



where Q!(x,t) = a{x,t) + ib(x,t). This is the leading order in the approximation: the correction to the 
amplitude b{x, t) of the solution q{x^ t, e) is purely a phase-shift in units of ^, where e^^ — the integer 
K counts the number of "ripples" as we explore the genus two region. So, finally, 

g(a;,i,e) = -2eT'f(°°) lim z {£{z)^ ^{z)),^ = 

Z—^OO 

= (^~b{x,t)e^'^'^ + 2{GK,^ + GK.>.)i2 + 0{V^)'^e^s{^^ K<x<K+l (B.30) 

This proves formula ( |7-1| . To verify the assertion regarding the different behaviors ( 7-3 ) we recall here 
equation (5-171 

GK,>t: + Gk.k = Fk+i,x: + Fk+i,>c + Ck + Ck ■ (B.31) 
We know that the matrices Fk+i,^^ Fk+i.^^ both are of order 0(e^^^^5) (by inspection of (5-11 1 and 



(5-13)). On the other hand 



Gk + Ck= lim z (*A'^'^.^+i) = 



-b{x,t) 



^i2K^ _ g-2(i<-+l)V 



(B.32) 



And this concludes the proof. Q.E.D. 
B.4 Proof of Thm. [721 

Looking at the formulae ( |5-10| and recalling formulae (4-271 and uk = — ma' we find 

-1 



A 



K 



:= (det[B'^,B^]B,,a'+det[B;^,B,]B,,2 -2 

|2||B 114 



2 ,||B ,|pB,,2B,„i 



Aa 



l + u, det[B;,B,,] 



K \ II A' I 



4(3(7,))^ 



> 0. 



(B.33) 
(B.34) 
(B.35) 
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At this point one needs to use ( A.3| (which defines the column vectors Ak,^k) and (A.7l (which gives 



A^,B'^), which produce (7-51. Plugging into the various object we obtain 

1 + IAoP, 



det[B^,BK] 



2\l 



2K 



B 



The rest is simply an algebraic manipulation which is as straightforward as uninteresting. Q.E.D 



\2K 



(B.36) 
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